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Abstract. We calculate the one- loop contributions to the MSSM effective potential when the 

scalar top fields have non-zero vacuum expectation values and study their impact on charge and 

CN ' colour breaking bounds. 

>■ 
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^3 I Supersymmetry is presently the best theoretical extension of the Standard Model (SM) of 

Cp ■ particle physics. Its popularity derives from the fact it solves both the hierarchy and naturalness 

Q . problems, stabilizing the weak scale by keeping huge quadratic divergences under control. Other 

"^ I incentives for supersymmetry include it making gauge unification more plausible and the tanta- 

O ' lizing prospect of supergravity theories aiding in the formulation of a quantum theory of gravity. 



The most unpleasant aspect of supersymmetry is the amount of free parameters it introduces: 
q"| assuming soft supersymmetry breaking we would have, in addition to the ~ 20 free parameters 

i-^ ' of the SM, about 100 new input values. Charge and/or colour breaking (CCB) bounds have 

been used widely to try to limit the size of that parameter space. Its foundations are simple: 

K> , in the Minimal Supersymmetric Standard Model (MSSM) electroweak gauge symmetry is bro- 

^ I ken when the neutral components of the two Higgs fields acquire non-zero vacuum expectation 

values. As these fields are also colourless the remaining gauge symmetry is SU{3)c x U{l)em, 
which is according to the experimental observation that we live in a world where electric charge 
and colour symmetries are unbroken. As we want to preserve angular momentum conservation, 
only fields of spin zero may have a non-zero vev. In the SM this is hardly a problem, since there 
is but one scalar field in the theory. 

In the MSSM we have a plethora of scalars, most of them carrying colour, charge or leptonic 
number. It is then theoretically possible that, should some of these fields have vevs, the MSSM 
potential develop a minimum deeper than the usual one, where only Hf and H2 have vevs. In 
that case the deeper minimum would surely violate charge and/or colour conservation. This 
appealing idea was first introduced in P and used to impose bounds on the supersymmetric 
input values: if a particular combination of them causes the appearance of a CCB minimum, 
we can exclude it on experimental grounds. A great deal of work was done on CCB bounds [2], 
of particular importance that of Gamberini et al. [H], showing the relevance of the one-loop 
corrections to the effective potential for gauge symmetry breaking, and the need to choose the 
appropriate renormalisation scale M to better take them into account. The choice of M became 
of fundamental importance in evaluating CCB bounds, so as to avoid an excessive constriction 
of the parameter space. A thorough review of the matter may be found in ref. 4_ and recent use 
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of CCB bounds in phenomenological studies in refs. [Hj- A common point in all of these analysis 
was that only the tree-level potential was studied. It was argued in ref. [SI that a careful choice 
of renormalisation scale would correctly reproduce the effects of the one-loop minimisation of 
the potential, both for the "real" minimum as for the CCB one. This led to comparing the value 
of the potential at the two minima at two different renormalisation scales. 

Partial one-loop contributions to the CCB bounds were considered in ref. ^ using renor- 
malisation group arguments. Recently a full one-loop charge breaking effective potential was 
calculated for the case where the scalar fields f^ and tr, not simply H^ and H2, acquire a 
non-zero vev [7]. At the same time it was argued that one should compare the values of the 
"real" potential and the CCB one at the same renormalisation scale due to the presence, in the 
effective potential, of field-independent terms. The minimisation of this one-loop CCB potential 
showed the Gamberini et al. strategy was not working for the CCB one-loop contributions [Hj 
and raised serious questions about the perturbative believability of bounds derived from this 
particular CCB direction. This would seem to be closely related to the smallness of the tau 
Yukawa coupling associated with this pattern of CCB, which generated typical CCB sparticle 
masses of the order of tens of TeV. CCB associated with a larger Yukawa should not fall victim 
to the same perturbative problems, and the only Yukawa coupling that satisfies those conditions 
is the top one ^. In this paper we undertake the calculation of the full one-loop effective po- 
tential for the CCB case where the scalar fields i^ and t/j, H^ and H2, acquire non-zero vevs. 
In section Q we review the model we will be working in and our conventions. We also discuss 
our strategy for CCB bounds and review the results of refs. [3IH1- The CCB mass matrices 
are introduced in section Q and the minimisation of the one-loop potential, and its results, is 
shown in section Q. We conclude with an analysis of the results and their relevance for CCB 
bounds. 



1 Charge and colour breaking in the MSSM 

We follow most of the conventions of ref. |9j (except for the sign of the // parameter). Because 
the Yukawa couplings of the first and second generations are so small compared to those of the 
third, we set them to zero. The MSSM superpotential is thus given by 

W = \tH2QtR + XbHiQbR + KHiLtr + ^iH2Hi , (1) 

with SU{2) doublets Hi = {H^ , Ff ), i^a = {H^ , H^), Q = {tL,hL),L = {vl , tl) and singlets 
iRi bR and tr. The MSSM tree-level effective potential is the sum of the F and D terms and 
the soft supersymmetry breaking potential. Its full expression may be found, for instance, in 
references |Zj and [H] . At the renormalisation scale M the one- loop contributions to the potential 
are 

where the Ma are the (tree-level) masses of each particle of spin Sa and n^ = (— l)2«a (2sa + 
1) Cq, Qa, Ca being the number of colour degrees of freedom of each particle and Qa counting 
its particle/anti-particle states. If only H^ and H2 acquire vevs vi/\/^, V2/V2, the negative 
contributions to the tree-level potential come from the — B ^viV2 term (and usually from the 
Tn^j V2 term as well). So that CCB occurs extra negative contributions are necessary - these 
come from terms cubic in the vevs, from the soft and F-term potentials. There are many 
possible choices of CCB directions, but one must remember that along with negative trilinear 
contributions to Vq come potentially large positive quadratic and quartic terms as well. In 

^The bottom Yukawa is larger than the tau one, but stiU of the same order of magnitude. 



refs. [ZllHl we studied in detail the case where the fields f^, tr had non-zero vevs. This particular 
direction, being associated with the tau Yukawa coupling, was in principle quite favourable to 
CCB - not only are rriL, rrir usually the smallest of the soft masses but the size of A,- should 
reduce the magnitude of the F-term contributions to the potential. Finally, as vevs at such a 
CCB minimum should be of the order of t; ~ 52 ^t/^t, the resulting potential, if negative, ought 
to be much deeper than the "real" minimum. The CCB direction we propose to study in this 
article - the scalar fields ti and i/j acquiring vevs - has none of these advantages: mq, mt are 
usually the largest of the soft masses; the top Yukawa being of order 1, the F-term contributions 
will be large; and the CCB potential, if negative, is not at all guaranteed to be deeper than the 
"real" one. Nonetheless, as explained in the introduction, we hope this CCB direction is not 
afflicted by the perturbative problems found in ref. ^. Also, the arguments against top Yukawa 
CCB rely on an intuitive analysis of the tree-level potential. Such analysis, however, is not 
possible in the case of the one-loop contributions, which are very complex. Specifying notation, 
we must remember that i^, Ir are in the 3, 3 representations of SU{3) respectively - each field 
has therefore three colour degrees of freedom. To limit the size of the SU{3) D-terms it is best 
to choose both vevs having the same SU{3) index - to simplify calculations we chose the third 
colour index ■^. We emphasize that this choice is not the most general CCB case associated with 
the top Yukawa, merely the one we expect will produce more interesting minima. So, let us 
consider that {i)^ , i^} have vevs {(7/\/2 , t/\/2} and {H^ , i/g} have vevs {fi/\/2 , V2/V^} as 
usual. The tree-level potential then becomes 

^0 = -f [viiq'^ + f) + q'^f] - -^{AtV2 - ^lVl)qt + -{mlvf + mlvi + mlq'^ 

+ fl (9^ - t'f , (3) 

with ml = m^ + fi^ and m2 = mj^ + //^. This potential being a 4- variable function 
analytical studies of CCB are possible only in simplified cases. The usual strategy Q^ considers 
only the tree-level value of the potential and the tree-level derived vevs: the latter simplification 
is based on ref. f^ where it was showed that for the "real" MSSM potential, by choosing a 
renormalisation scale M of the order of the largest mass present in AT/i, the tree- level vevs were 
a good approximation to the one-loop ones - that is to say, in this range of M the one-loop 
contributions to the minimisation conditions, given by 

^ 327r2 ° dvi V M^ J 

are not significative. Notice, however, that this does not necessarily imply that the contributions 
AVi will be negligible. In fact, even in the usual MSSM calculations with an appropriate mass 
scale, the value of the one-loop potential is positive whereas the tree-level one is negative. This 
is not a problem since the value of the potential, CCB bounds excluded, has no bearing on the 
phenomenological aspects of the model. But it shows that even if an adequate choice of scale 
may render insignificant the derivatives of the one-loop contributions to the potential, that same 
choice does not imply that AVi itself is negligible. Since the typical mass is usually different in 
the "real" and CCB potentials this leads to both being compared at different renormalisation 
scales. It was this point that led us to argue against this procedure in ref. [7], based on the fact 
that the sum Vq + AVi is not renormalisation group (RG) invariant. Instead, the complete RG 

^Identical results would be obtained for the colour indices (1) or (2). 



invariant effective potential is given by jEIl 

V{M,\i,4>j) = ^{M,Xi) + VoiXiAj) + hAVi{M,XiAj) + o{h^) , (5) 

wliere Aj stands for the couplings and masses of the theory and (t)j for its fields. The field- 
independent function Jl, implicitly or explicitly, depends on the renormalisation scale M. The 
only difference between the CCB and "real" potentials being the different set of values some 
of the fields 4)j have, Q. is the same in both cases - so, if we compare V^^^^^'^ and V^'^^ 
at different scales the contributions from Q will not be considered correctly. But comparing 
potentials for the same value of M means that at least in one case the one-loop contributions 
to the vevs will have to be included. By consistency, if h contributions to the vevs are being 
taken, h contributions to V, that is, AVi, must also be considered. In ref. jS] we found an even 
stronger argument for studying one-loop CCB potentials and their minimisation: eq. © implies 
that d{Vo + /\Vif^^/dM = d(Vb + IWi)^'^^^^^ /dM, up to two-loop effects. In other words, 
the two potentials must run with the renormalisation scale parallel to one another. In we 
found that the CCB vevs were so large that the one-loop CCB effective potential was not RG 
invariant. Perturbation theory had broken down and the typically small two- loop effects had 
become quite large. As a consequence, the condition "y^'C'^ < yMSSMn ^^^ renormalisation 
scale dependent and as such bounds thereof derived not reliable. We expect this will not happen 
for the top-associated CCB direction, as the typical vevs should be smaller. 

In short, our CCB strategy will be to calculate the one-loop CCB potential, obtain from its 
minimisation the one-loop vevs {vi , U2 , <? , t} and compare it to the one-loop minimised "real" 
one-loop potential, at the same renormalisation scale. For comparison, we will also perform 
tree-level minimisations of the potentials (both MSSM and CCB) and compare their values - 
even though, as follows from the discussion above, this tree-level procedure is misleading and 
may induce errors. For the one-loop calculations we will need the CCB masses contributing to 
AVi. The complication that arises, as in ref. [Jj, is the fact that the vevs q and t cause the 
mixing between charged/neutral and coloured/colourless fields in the theory. For instance, the 
trilinear terms in the soft potential cause mixing between {H^ , H2 , i^ , i^ } - and further 
mixing between them arises from the F and D-terms. Likewise, similar mixing occurs between 
{H^ , H2 , i^ , ins- The charged, pseudoscalar and CP-even Higgses will therefore have 4x4 
mass matrices. Because colour symmetry has been broken, particles carrying colour indices 
{(1) , (2)} will have different masses from those with colour index (3) - the {(1) , (2)} squarks 
have mass matrices very similar to those of the non-CCB case. The existence of vevs carrying 
colour degrees of freedom gives mass and electric charge to four gluons, three others remaining 
massless. The eighth gluon remains neutral but becomes also massive by mixing with the -B^ 
and WJ^ fields. For the fermions the mixing is even more extensive: the charginos become a 
mixture of charged SU{2) gauginos, the fermionic partners of the charged Higgses and the (3) 
component of the bottom quark. The case of the neutralinos is even more complex: a 7 x 7 mass 
matrix originates from the mix of neutral U{1) and SU{2) gauginos, fermionic partners of the 
neutral Higgses, (3) components of the top quark and the eighth gluino. The gauge interaction 
term between quarks, scalar quarks and gluinos, —igs (jr Aj 4)G^ /\/2, also causes mixing between 
the {(1) , (2)} quark components and the G^''^ gluinos - the results are two 4x4 identical mass 
matrices. The mass of the gluinos G^-'^ remains unchanged, M3. We present the CCB mass 
matrices in the next section. 



2 Mass matrices 

Let us define the coefficients G as 

r-i 2 2,^2 4j.2 ri 2 2 2 ni 2 j.2 i a\ 

Gi = ^2 - wi + -g - -t , G2 = V2 - v^ - q , G3 = q - t . (6) 

We now list the masses of the MSSM when q and t are non-zero. 

• First and second generation sleptons and sneutrinos (rii = n2 = Ui,^ =2x2) ^: 

Ml=m% - ^—Gi + |g2 Ml=ml + ^—G^ 

Ml=ml-'-^G,-'^G, . (7) 

• First and second generation squarks, colour indices {1,2} (ni = n2 = 2 x 4, for both up 
and down type squarks): 

A^2 (1:2) _2, 9'^ ^ 9I ^ 9I ^, , .2 (1:2) _ ^2 5'^ ^ I ^3 ^ 

7W^j - m^ + — Gi - — G2 - — G3 7W^2 - "^u 6~ 12 

M|(''')=m|> + ^Gi + ^G2 - ^Ga mJ ^^''^ = m^ + ^ Gi + ^ G3 . (8) 
di ^ 24 8 12 '^a '^ 12 12 ^ ^ 

• First and second generation squarks, colour index 3 {n\ = n2 = 2 x 2, for both up and 
down type squarks): 

M|/^)=m| + |^Gi-|G2 + |G3 M, 

M|/'^=m| + ^Gi + |G2 + |G3 M|;-=m^ + ^Gi-^G3 . (9) 

• Tau sneutrino and sleptons {rii,^ = rii = 722 = 2): 

Ml=ml- ^Gi - |g2 , {Ml\ = h ^] , (10) 



r2 (3) 

M2 


2 9'' 

- m„ g <^i - 


9I ^ 


r2 (3) 

d2 


= ml + ^G^- 


93 r< 



with 



)2 /2 2 ^ 

Of = mi + ^ uf - — Gi + -|- G2 6f = -| (/ii;2 - ^r I'l) 

Cf = ml + ^vl + ^-^G^ . (11) 

• Charged (n = 6) and neutral (ni = n2 = 3) electroweak gauge bosons: 

''"'GO bco dco 

boo CqO 6^0 

ydco eQO fco 



1 
4 



Mw = 1 9I {vl +vl + q^) , [Ml,] = ( boo Cgo 6^0 | , (12) 



with 



52/2, 2, 2\ u 92 9 I 2 , 2 ^ 
Ago = ^(^^1 +V2+q ) 000 = — {V^+V2- -( 



We considered degenerate first and second generation sparticles, but these results are trivially generalised. 
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' / 1 



fco 



92 93 2 



(13) 



This matrix has one zero eigenvalue, corresponding to a "photon" resulting from gauge symmetry 
breaking ^. 



Charged gluons (ni = 4 x 3): 



MS. = ^5i(9^ + t^) 



(14) 



Top scalars, colour indices {1,2} (ni = 722 = 4): 



[^f 



2(1,2)1 



with 



1 x2„.2 , 5' 



,/2 






ffi 



(15) 



9i 



(^i = m'Q + -A^2^ + ^Gi - ^G2 + f^(3(z^ - G3) 



24 



12 



c- = m? + i A? «i - ^ Gi + ^(3t2 + G3) . 
t t I 2 * 2 6 12^ ^ 



(16) 



• Bottom scalars, colour indices {1 , 2} (ni = n2 = 4): 

r2(l,2)i 



[M 



% h 



h % 



(17) 



with 



1 o'2 ^2 2 

2,^\22,y,^ , 92 ri 93 r^ 

ai = ^Q + 2^bVi + ^ '^1 + -§- ^ ~ 12^ 



^6 



-j={AbVi - HV2) 



1 x2„,2 , 9' 



9l 



Cfo = "ife + o A^ vf + ^ Gi + ^Ga . 



12 



12 



(18) 



Charged Higgs (mix between H^ , H2 , b\^ and b^ , ?^l...4 = 2): 



/a± b± d± e±\ 

b± c± f± g± 

d± f± h± i± 

\e± g± i± j±J 



(19) 



with 



1 /2 2 

2,-'-\22 9 ^ , 92 / 2 , 2 2\ 



^The symmetry breaking we have chosen leaves intact a SU{2) x t/(l) gauge group, corresponding to an integer 
charge quark theory. In this theory four gluons couple directly to the photon and as such possess electric charge. 



See reference |11| for details. 



b± = B fi + ^viV2 

1 /2 2 

c±=m2 + -Xtt + —Gi + —{v^+v-^+q) 
d± = y=nt + -{gl-2Xl)viq 

e± = -\bXtV2t T^'^bQ 

f^ = ^Att+^{gl-2Xf)v2q 

9± = --T^fJ'Q + ^XbXtVit 



./2 



92 /„,2 



9l 



i± = —={Ay,Vi - IIV2) 

1 '2 

2|^\2/2| 2\ , 9 ^ 

j±=mf, + -Xf,{v^+q) + -^Gi 



9l 



G, 



(20) 



tJ-O +(3) 



r(3). 



Pseudo scalars (mix between the imaginary parts of H^, H2, t^ and tj^ ; ni...4 = 1): 



[M 



^qJ 



^H ^H fu 9h 

dn fn hfj ijj 



(21) 



with 



2 9 ^ 92 ri 



as = ml + '-X^g' + t') + '-^G, + 'Ig2 



9h 



^/2 
^/2 



Atq 



hu = ^% + ^A?(.2^ + t2) + Cgi - #G2 + f G3 



/2 

24 



53 
6 



V2 



ijj = ^=(AtU2 - iivx) 



1 /2 2 

2|-'^\2/2i 2\ 9 ri 9'i ri 

3H = ^t + 2 ^t (^2 + 9 ) - -^ Gi - -^ G3 



(22) 



[7-0 l(3) 



Higgs scalars (mix between the real parts of H^, H2, t)^ and r^ , ni...4 = 1): 



7(3). 



[M^o] 



bn CH fn QH 

dn Jh hn in 

\eH QH in JH/ 



(23) 



with 



aH = ml- ^{Gi-2vl) - ^{G2-2vl) 

o o 

bn = -B fi - -{g +5-2) n V2 

1 /2 2 

CH = mi + - A? (g^ + *') + ^ (Ci + 2t;i) + ^ (G2 + 2t;i) 

At ^ 1 /2 , 
en = ^^9+25 vit 

fH = >HV2q j=Att + —{g - 85-2) V2 q 

1 



5if = - — I ^t 9 + ^ (3 At - g'^) V2 t 

hH =ml + IxUvl + t') + '^ (g, + Iq' 



f (G2-2g2) + |(G3 + 2g2) 



^H 



^(AtV2 - ^.v^) + ^{9X', - g" -3gl)qt 



1 



2 I -^ \2 / 2 I 2\ 

jH = mt + -Xt{v2+q ) 



^(G,-?*^)-fG,-2.^, 



(24) 



Tau lepton and bottom quark, colour indices {1,2} (n^ 

\rVi/V2,Mb = XbVi/y/2. 



-4, Ub = -2 X 4), M^ 



Charginos (mix between the fermionic partners of the charged Higgses H^ , H2 , the SU(2) 

b^^^ and 6^^; 



gauginos ^~^ , ^~ and the colour index 3 components of the bottom quark, b\ and br 



ni..s 



-2 : 



[M^±] 
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M2 



-^-1 



V 








-M 







t 



M2 

ft -2 






g2 

V2 



■ui 









V2^ 






V2^ v^ 



I'l 



\ 







^ 



vi 



J 



(25) 



• Neutralinos (mix between the fermionic partners of the neutral Higgses H^, H2, the U(l) 
and SU(2) gauginos B, VF3, the colour index 3 components of the top quark, t\ , tj^ and 



eighth gluino G 


; ?^l...7 = 


-2): 














( Ml 





-^vi 


^V2 


-i^ 


29'^ 
3 *- 










M2 


f^i 


-fv2 


fq 










-^vi 


fvi 





-fj- 











[M^o] = 


f- V2 


-fv2 


-M 





V2* 


-7l^ 







-i^ 


fq 





-** 





V-2^^ 


.93 




'-ft 








-al' 


^^2 





.93 f 




[ 











-^'^ 


^* 


Ms 



(26) 



Gluinos (mix between the colour index {1,2} components of the top quark and the gluinos 



(74-7. „^ ^ = -2 X 2 5) 



[M, 



Gi 



( A^2 -^9 \ 



V 



%V2 

-^. 



^* 
Ms 

Mg 



(27) 



The remaining gluinos (G^''^) contribute to the one-loop potential with mass M3 and are affected 
by a factor n = —2 each. At the tree-level minimum the matrices p5|) and ()19p have each a zero 
eigenvalue, the matrix (|21|) has two. Counting the multiplicities of each particle this corresponds 
to a total of eight Goldstone bosons, corresponding to eight gauge bosons that acquire mass (the 
Z*^, the W and five of the gluons). We checked these mass matrices by computing Str M'^ = 
J2a ™a M^ - because supersymmetry is broken in a "soft" way this quantity should be field- 
independent (that is, it should be independent of the value of the vevs), and it is simple to 
verify that condition is met. 



3 Minimisation of the one-loop CCB potential 



With the mass particles for the whole sparticle spectrum obtained it is a simple task to obtain 
the derivatives of eq. ^ and perform the one-loop minimisation of the CCB potential as was 
done in ref. jH|. Here we will use a different method to minimise the potential: we will not bother 
with its derivatives and instead look directly at the values of Vq + AVi as function of {vi ,V2,q, t) 
by means of a modified simulated annealing algorithm incorporating Q-sampling of the phase 
space |13j . The code for the application of this algorithm was developed by J.M. Pacheco. The 
reason for choosing this method is purely a practical one: the algorithm is extremely efficient and 
the computation time necessary for a scan of the MSSM parameter space drastically reduced. 
We checked the code by comparing its results with those of a MSSM calculation using the 
potential's derivatives, for a vast region of parameters. We also compared its results with those 
obtained using the numerical minimisation tools of the MATLAB package. All these approaches 
produced the same results. 



5.(1) .(1) 



fj^' mix with G 



{& + iG"')/%/2, G" = (G^ - iG'*)/^/2. An identical mixing occurs between t'-'^'> , G^ 



and G producing degenerate gluino masses - thus, the extra factor of 2 in the coefficients n. 



In this work we will mostly study the MSSM with universal input soft parameters. At 
the gauge unification scale we choose common values Aq, Mq and mo for the A parameters, 
the soft gaugino and scalar masses. We took —4 < Aq < 4 TeV and 20 < Memo < 900 
GeV. Further, we have taken 2.5 < tan/3 < 10.5 and considered both possible signs for the fi 
parameter. This selection of parameter space is by no means an exhaustive one (we left out 
high values of tan/3, for instance) but is already very extensive and should provide a good 
idea of the importance of the CCB bounds derived from our one-loop potential. We follow the 
"top-bottom" approach outlined in ref. [H]: at the weak scale Mz we input Mz = 91.19 GeV, 
Mt = 167.2 GeV, Mf, = 2.95 GeV, M^ = 1.75 GeV (these are running fermion masses, not the 
pole ones), ai = 0.01667, 02 = 0.032, as = 0.1 and the value of tan/3. Because the gauge and 
Yukawa /3- functions do not depend on the soft parameters ^, we can run these parameters up 
in the energy scale until we find the gauge unification scale, defined as the value Mjj for which 
the couplings ai and 02 meet. At that point we input the values of the soft parameters, chosen 
as explained above, and run the whole theory to a scale Mq = max(M^ , Mq , 53 Ac/Xt) - this 
scale is a good estimate of the heaviest masses in the theory, thereby reducing, in principle, the 
size of the logarithmic contributions to the potential. At that scale we use the one-loop MSSM 
minimisation conditions (see, for instance, ref. 0) and determine - if possible - the parameters 
/i and B. We then calculate the sparticles' masses ^ and use recent experimental bounds ^^ 
to reject those "points" already in contradiction with observational evidence. We end up with 
over 39000 "points" of parameter space. We also studied a small non-universal parameter space 
where we took the universal "points" obtained earlier and set the m^ soft parameter to negative 
values (of the order of, at the most, (100 GeV)^), as this situation would seem one of the likeliest 
to result in CCB. This situation clearly requires non-universality, as universal values of the soft 
masses are unlikely to result in negative m^ at the weak scale, due to the form of the /3-function 
for this soft parameter. With the new value of m^ we minimised the (new) potential and, again, 
used experimental sparticles' mass bounds to reject "points" in disagreement with observations. 
Our non-universal parameter space ended up with about 4000 "points" . 

We set out to determine the impact of the one-loop contributions on CCB bounds and as 
such we endeavored to compare results coming from tree-level and one-loop minimisations of the 
potential. The tree-level minimisation of the MSSM potential is performed analytically in the 
standard way, see for example [H]. So, both at tree-level and one- loop, we compute the value 
of the MSSM potential, V , and perform the numerical minimisation of the CCB one, 

obtaining the vevs {vi , V2 , q , t}. With these we calculate the value of y^'^^ and compare it 
with y'^''^SS'^'\ Finally, a word on thresholds: we follow the procedure of ref. |^, using the full 
MSSM /?- functions from Mz to Mjj, and choosing the input parameters at Mz such that the 
threshold contributions are automatically taken into account - this is an effective procedure, but 
shown to produce good results. And for our purposes - determining if CCB occurs or not - this 
degree of precision should be more than adequate. Being based on a Monte Carlo method the 
algorithm depends on the initial conditions used, so several runs of the program were necessary. 
We only accepted those extrema with all CCB squared masses positive, except the expected four 
zero eigenvalues corresponding to the Goldstone bosons. We remember that when performing 
a one-loop minimisation of the MSSM potential it is usual to find negative squared masses in 
the Higgs sector [2j. They correspond to the Goldstone bosons which have zero masses for a 
tree-level minimisation. Since we do a one-loop minimisation and compute their masses using 
tree-level matrices, negatives do occur sometimes. However, the absolute value of these negative 
squared masses is very small when compared to the other masses in the theory and thus, based 
on ^H], they can be safely set to zero. 

^Except indirectly, in the form of particle thresholds. 

^With the tree-level mass matrices except for the neutral CP-even Higgses, for which full one-loop expressions 
were used. 



The results of this scan of the MSSM parameter space may be resumed as follows: 

• No unbound-froni-below (UFB) directions were found for the one-loop minimisations. 
They appear frequently, however, for the tree-level minimisations. It is easy to see from 
equation (jSJ, for instance, that with q and t constant and ?;i ~ W2 — ^ oo we have a potential 
tree-level UFB direction. These directions are characterized by vevs assuming arbitrarily 
large values, causing the potential to become arbitrarily negative. No such thing happened 
in the one-loop minimisation of the potential, which confirms the expectations of the au- 
thors of s4j : they had assumed that the UFB directions they had found at tree- level might 
not be present if a one-loop minimisation was performed. We confirmed that that is the 
case, and that the one-loop contributions to the potential stabilize the vevs (which is also 
in agreement with the conclusions of p!). 

• For the tree-level minimisations, the value of the CCB potential was found to be always 
deeper than the MSSM one. For many of these minima, however, the values of q and 
t found were practically zero. What happens is that the tree-level MSSM potential is 
found for very specific values of the vevs (namely such that vf + V2 = (246 GeV)^ and 
tan/3 = ^2/^^!)) but that is not the only possible MSSM minimum - the tree-level MSSM 
minimisation conditions consist of two coupled cubic equations and as such can have as 
much as nine different solutions. What we find is thati, for the CCB case, we let the four 
vevs "roam freely" as we minimise numerically the potential and, as such, many solutions 
with g ~ t ~ are found. However, those solutions have values of f 1 and V2 which are 
not phenomenologically acceptable (giving as they do wrong values for the gauge bosons' 
masses, for instance). These alternative minima, which preserve the MSSM gauge sym- 
metry, are actually deeper than the "standard" minimum. But again we must remember 
that a comparison of tree-level potentials with tree-level derived vevs is inherently flawed, 
and we only undertook it to show the radical differences with the one-loop case. 

• The MSSM potential is always deeper than the CCB one at one-loop, both for the universal 
and non-universal parameter spaces. There are now no "alternative MSSM minima" as 
in the tree-level case, and the reason is easy to understand: while a given combination 
of supersymmetric parameters might correspond to several possible combinations of vevs 
as minimisation solutions at tree-level, it is almost impossible that that happens for the 
one-loop potential, given its complexity. 

The comparison between these two final points reveals how different a one-loop procedure can 
be from the tree-level one. But it also requires some explanation: what is it about the one-loop 
CCB contributions that raises the value of the potential, always above its MSSM value? First, 
we have already seen the importance of a one-loop minimisation to stabilise the values of the 
vevs (avoiding UFB directions and preventing the appearance of "alternative MSSM minima"). 
Secondly, the consequences of the different gauge symmetry breakings become apparent only 
at the level of the one-loop potential, when we consider the different spectra of masses thereof 
resulting. In section |21 we showed that a very significant difference between the MSSM and the 
CCB case is the mass of the gluons: all zero for the MSSM, five of them gaining mass for CCB. 
For CCB vevs q and t of the same order of the electroweak vevs vi and V2 (as we expect them 
to be, if not even higher), we expect the gluons to have masses considerably larger than Mz 
or Mw since, from eq. dJ, their masses are proportional to the strong coupling constant. As 
these are bosonic masses, their contributions to AVi, from eq. ((2)), should be positive. Assuming 
a CCB minimum exists, then, its gluon one-loop contribution to the potential is expected to 
be large and positive, whereas for the same SUSY parameters the similar contribution to the 
MSSM potential is zero. This is a possible explanation, but the one-loop contributions are very 
complex and the real reason may be lying with other terms. 



Ultimately, all we can conclude from these results is that we have a "numerical demonstra- 
tion", for this vast parameter space, of the MSSM minimum being deeper than any possible CCB 
minimum. This follows similar conclusions reached in ref. |^, where no acceptable CCB minima 
associated with the tau Yukawa coupling were found. Also, recent work ^3] demonstrated that 
for two Higgs doublet models (2HDM) at tree-level, if a minimum preserving U{l)em and CP 
symmetries exists, it is a global minimum. In such models it is impossible to have tunneling for 
charge or CP breaking vacua. Of course, in this work we are dealing with charge and colour 
breaking, in a model that is far more complex than the 2HDM. We are also studying one-loop 
minima, not tree-level ones. But these different elements may be suggesting that the occurrence 
of CCB is not as easy to occur as a tree-level analysis in the MSSM leads us to believe. 

It is important to remark that we did not exhaust the MSMM parameter space. The pos- 
sibility of dangerous CCB minima occurring may happen in some portion of parameter space 
not included in our study. Even though the range and number of parameters we chose was very 
general, this possibility cannot be wholly dismissed. Finally, all our calculations have assumed 
that the MSSM vacuum is the absolute minimum of the theory. There is however j2j the the- 
oretical possibility that the the CCB minimum is deeper than the real one and the tunneling 
time between both vacua be superior to the age of the universe. This leads to a relaxation of 
the CCB bounds usually obtained. In our case we do not need to worry about this possibility 
since no CCB absolute minima were found. 

We hope to have convinced the reader of the importance of full one-loop calculations in 
estimating CCB bounds. In references [7j, |Sj and this work we showed that for RG consistency 
one had to compare the one-loop MSSM and CCB potentials, and that the results of that 
enterprise gave very different results from tree-level studies. We studied very specific CCB 
directions and vast sections of the MSSM parameter space, and no CCB minima deeper than the 
"normal" vacuum were found. Although at this stage we cannot make sweeping generalisations 
and claim all CCB bounds used in the literature are wrong, these results urge some caution. 
In reference ^, for example, large sections of parameter space are excluded on CCB grounds, 
and predictions of supersymmetric masses are affected by it. If indeed, as this work suggests, 
tree-level CCB bounds are over-estimated, we may be excluding areas of parameter space that 
can be of experimental interest. A careful re-evaluation of CCB bounds might be in order. 

Acknowledgments: My deepest thanks to Jorge Pacheco for providing me with the code 
for the minimisation algorithm. This work was supported by a fellowship from Fundagao para 
a Ciencia e Tecnologia, SFRH/BPD/5575/2001. 



References. 

[1] J.M. Frere, D.R.T. Jones and S. Raby, Nud. Phys. B222 (1983) 11. 

[2] L. Alvarez-Gaume, J. Polchinski and M. Wise, Nud. Phys. B221 (1983) 495; 
J.P. Derendinger and C.A. Savoy Nud. Phys. B237 307; 

C. Kounnas, A.B. Lahanas, D.V. Nanopoulos and M. Quiros, Nud. Phys. B236 (1984) 438; 
M. Claudson, L.J. Hall and I. Hinchliffe, Nud. Phys. B228 (1983) 501; 
M. Drees, M. Gliick and K. Grassie, Phys. Lett. B157 (1985) 164; 
J.F. Gunion, H.E. Haber and M. Sher, Nud. Phys. B331 (1988) 320. 

[3] G. Gamberini, G. Ridolfi and F. Zwirner, Nud. Phys. B331 (1990) 331. 



[4] J.A. Casas, A. Lleyda and C. Munoz, Nud. Phys. B471 (1996) 3. 

[5] U. Ellwanger and C. Hugonie, |hep-ph/9811386" 

S. Abel and T. Falk, Phys. Lett. B444 (1998) 427; 

S. Abel and C. Savoy, Phys. Lett. B444 (1998) 119; 

S. Abel and B. Allanach, Phys. Lett. B431 (1998) 339. 

OPAL Collaboration, Eur. Phys. Jour C7 (1999) 407; ibid., Eur. Phys. Jour C12 (2000) 
567. 

[6] H. Baer, M. Brhlik and D. Castano, Phys. Rev. D54 (1996) 6944. 

[7] P.M. Ferreira, Phys. Lett. B509 (2001) 120; Err. Phys. Lett. B518 (2001) 333. 

[8] P.M. Ferreira, Phys. Lett. B512 (2001) 379; Err. Phys. Lett. B518 (2001) 334. 

[9] V. Barger, M.S. Berger and P. Ohmann, Phys. Rev. D49 (1994) 4908. 

[10] C. Ford, D.R.T. Jones, P.W. Stephenson and M.B. Einhorn, Nucl. Phys. B395 (1993) 17. 

[11] P.M. Ferreira, hep-ph/0210024i 

[12] J. Rosiek, Phys. Rev. D41 (1990) 3464. 

[13] A.K. Hartmann and H. Rieger, Optimization algorithms in physics, Wiley VCH, Berlin 
2002. 

[14] K. Hagiwara et al, Phys. Rev. D66 (2002) 010001. 

[15] P.M. Ferreira, I. Jack and D.R.T. Jones, Phys. Lett. B387 (1996) 80. 

[16] D. Pierce, hep-ph/9407202| 

[17] P.H Chankowski, Z. Pluciennik and S. Pokorski, Nucl. Phys. B439 (1995) 23. 

[18] Y. Fujimoto, L. O'Raifeartaigh and G. Parravicini, Nucl. Phys. B212 (1983) 268. 
E.J. Weinberg and A. Wu, Phys. Rev. D36 (1987) 2474. 

[19] P.M. Ferreira, R. Santos and A. Barroso, Phys. Lett. B603 (2004) 219. 



